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Abstract. The Rosette Complex in the constellation of Monoceros is a magnificent 
laboratory for the study of star formation. The region presents an interesting scenario, 
in which an expanding HII region generated by the large OB association NGC 2244 is 
interacting with a giant molecular cloud. Inside the cloud a number of stellar clusters 
have formed recently. In this chapter we present a review of past and present research 
on the region, and discuss investigations relevant to the physics of the nebula and the 
molecular cloud. We also review recent work on the younger embedded clusters and 
individual nebulous objects located across this important star forming region. 



1. Historical Perspective 

The Rosette Complex (/ = 207.0, b = —2.1) is located near the anti-center of the 
Galactic Disk in the constellation of Monoceros. The region is very popular, partly 
because of the staggering beauty of its main feature: a very extended emission nebula 
which hosts a large central HII region, evacuated by the winds of a central OB associa- 
tion (see Figure [B. 

The complex is part of a much larger structure known as the Northern Monoceros 
Region. This region comprises the Mon OB 1 Cloud (host of NGC 2264 and the Cone 
Nebula), the Monoceros Loop (NGC 2252), and the Mon 0B2 Cloud in which the 
Rosette is one of the most prominent features (see Figure |2ll. 

The catalog name for the Rosette can be somewhat confusing because it is not 
unique: The nebula itself is usually cataloged as NGC 2237 or NGC 2246 (especially 
by amateur observers), although NGC 2237 originally referred to the brightest patch at 
its western side and NGC 2246 originally pointed to a bright zone at the eastern side. 
In addition, while the central cluster is usually known as NGC 2244, it has also been 
cataloged as NGC 2239. However, this latter designation historically referred to the 
brightest star in the region, 12 Monocerotis. 

The cluster was first noticed by Flamsteed in the late 17th century and later re- 
ported by William Herschel — who did not notice the nebulosity — and John Her- 
schel, who disco vered several co nspicuous nebulosity features and reported them in his 
general catalog (lHerscheliri864 NGC 2239 = GC 1420). A contemporary report on 



nebulous emission was done by Albert Marth (lLasselllll867l . NGC 2238 = GC 5361 = 
Marth 99), who described a "small, faint star in nebulosity." 
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Other parts of the nebula (NGC 2237 and NGC 2246) were reported by ISwift 



(|1885|) who cataloged the object as being "pretty bright [pB], very, very large [vvL] 
and diffuse [diff]." Afterwards, the region was formally known as the "Swift Nebula," 
until the name "Rosette" became more popular. The total exten t of the Rosette was not 
determined until the first photographic plates were obtained by Barnard' (1894). 

Two of the first applications of Rosette Nebula photographic data were made 
bv iHubbld (Il922h in his study of diffuse nebulae associated with massive stars, and 
Minkowskil (11949'). who published a photographic study along with a first discussion on 
the expansion of the HII region by the winds of O stars in the central cluster. Minkowski 
also pointed out the existence of elephant trunks and dark globules. He estimated the 
mass of the nebula to be W^Mq and suggested that it could be "surrounded and pro- 
bably embedded in obscuring material," thus proposing the existence of the companion 
molecular cloud. 

The first formal r adio wavelength observation of the Rosette Nebula was reported 
by lKo & KrausI ( 1955b . They surveyed the area with the Ohio State University antenna 
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Figure 2. The location of the Rosette in the context of the Monoceros Complex, 
from Perez (1991) 



at 187.5 MHz with a HPBW resolution of 1.3', and calculated a blackbody temperature 
of 200 K for the nebula. Later, Menon (1962) used the 85 ft. Tatel telescope to make a 
very complete study of the nebula at 3 GHz (10 cm). His analysis showed convincingly 
that the Rosette Nebula is an ionization bounded Stromgren sphere and he discussed 
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how it could have formed from a fraction of a large dark cloud with a total mass of 
approximately 1.4x10^ Mq. 



2. The Rosette Nebula and the Young Cluster NGC 2244 



2.1. The Rosette Nebula 



Qgura & Ishidal ( 1981 ) tabulated the different methods used to determi ne the age of the 
nebula. The se varied from studies of the properties of the central cavity ( Kahn & Menon| 



19611 : iLaskeri ,1966.) to evolu tionarv models of the HII region based on the luminosity 



of the star s faielhning| |l968h. Other methods involve time scales of radiation pressure 
(Mathews 19661. ~967h . estimates of the formation time for dark globules in the nebula 
(IHerbig 1974). and the separation of [OIII] emission lines (.Smith. 1973) . The median 
value of all these age estima tes is approximately 3 ± 1 x 10^ yr. 

A series of studies by ICelnikI (Il983l 1 19851 . 11986^ discussed the global physical 
characteristics of the Rosette Complex. The first two of these are dedicated to the 
nebula, while the third one is a model of the interaction with the molecular cloud, thus 
we defer its discussion to Section 3. In the first of the ICelnikI articles, he presented 
a map of the Ha emission in the nebula region, and calculated a total integrated flux 
density of 5 x 10^^^ W-m^^ within 60' from the center of the HII cavity. He suggested 
that the emission is contained in a more or less symmetric ring with a peak at 16' from 
the center. 

Table 1. Distance Estimates to the Rosette (NGC 2244) 



Author 


Value (pc) 


Method 


Johnson (1962) 


1660 


Photoelectric Photometry 


Ogura&Ishida(1981) 


1420 


Visual Photometry 


Perez et al. (1987) 


1670 


Visual Photometry 


Park & Sung (2002) 


1660 


Visual Photometry 


Hensberge et al. (2002) 


1390 


Spectroscopy 



In the second paper. ICelnikI reported radio continuum observations (1410 and 4750 
MHz) from which he was able to determine that the nebula is bound by ionization, 
forming a spherical shell with radius of ^^20-30 pc (in agreement with the 3.7, 7.7 
and 14.3 GHz observations of .Kaidanovskii , 1 980.) and a tot al ionized ma tter mass of 
2.3 X 10"^ M0 (almost twice the estimate of 1.1 x 10"^ M© by lMenon|[l962l) . Using the 



HI 12a and Hell2a (4619 and 4621 Mhz) radio recombination hues (RRL) Celnik 
calculated a He+ abundance of 0. 12ib0.03 and a non-LTE average electron temperature 
for the nebula of Te = 5800 ± 700 K — almost 1100 K above the LTE — with no 
evidence for a significant gradient with respect to the radial distance from the center. 

However, other studies reve aled that the LTE temperatu re of the Rosette was def- 
initely not uniform. For example iPedlar & Matthewsl (Il973 ) found Tg values between 
6400±1300 and 8160±1100 K (ave rage of 790 K) across the nebula using H166a 
RRL (1424 MHz). In a similar study. IViner et all fl979 ) calculated non-LTE values be- 
tween 3900±700 and 6000±500 K from HlOOoRRL (6480 MHz) observations. The 
25 MHz observat ions (decameter l ength) of iKrymMnl (Il978h suggested a median Tg 
of 3600 K. Later, iDeshpande et al.l (119841) using continuum absorption observations at 
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34.5 MHz combined with 2700 MHz observations by lGraham et alj (1982) argued that 
the electron temperature might vary from approximately 5000 K in the southeastern 
quadrant of the nebula — which contains more dust from the interface w ith the molec- 



ular c loud — to roughly 8000 K at the northwestern regions. Finally, iTsivilev et al. 



(l2002h presented observations of the H92q! RRL at three points coincident with crucial 
positions of the 4750 MHz map of Celnik which revealed a higher average tempera- 
ture of 7980±580 K. They also discussed all other RRL measurements of the Rosette 
Nebula available in the literature, and found Tg measurements to vary from 3000 to 

10050 K. 

The investigation by IShipman & Clc^ (1994) revealed a good fit to a T oc r^°^, 
a = 0.4 model for the temperature gradient in the nebula cavity which, interestingly, 
could not be adjusted to the observed IRAS emission. Instead, they found that this 
temperature gradient was better adjusted to a = 0.05 for r < 47' and a = 0.2 for 

47' < r < 65' (see also Section 3.2)^ 

The pioneering observations byl Gosachinskii & Khersonskiil ( 1982 ) of the HI emis- 
sion at 2 1 cm made with the RATAN-600 telescope revealed that the Rosette Complex 
and the Monoceros Loop were enclosed in a thin HI envelope 130 pc in diameter, possi- 
bly expanding at 20 km-s~^, suggesting that the supernova event t hat gave or igin to the 
Monoceros remnant occurred about 1.8x10^ yr ago. The study of iGuseva et al. (1984) 
discussed in detail how the supernova event in Monoceros could in fact be responsi- 
ble for the collision of clouds that triggered star formation in NGC 2264 (1 Mon) and 
the Rosette Complex (2 Mon). Later, Kuchar & Bania ( 1993) made a new, local 21 cm 
map of the Rosette Complex using the Arecibo telescope. They found that neutral gas in 
the Rosette Complex is distributed in three main regions which form a rough, extended 
shell of about 45 pc in radius around the nebula and extend beyond the molecular cloud. 
This shell would have a center of expansion at (a, 5) = (6''31™48^ 4'='59'12", J2000), 
and has a mass close to 2x 10^ M©, which implies a budget of kinetic energy for the 
shell expansion of approx. 4x10^^ ergs, or roughly 2% of the total energy available 
from the stars in NGC 2244. 



2.2. NGC 2244 

The prominent OB association NGC 2244 contains more than 70 high mass sources 
with spectral types O and B (7 and 24 sources respectively, see Table |2ll. The stars in 
the cluster are pre sumed responsible for the evacuation of the central part of the neb- 
ula. iFlvnnI (Il965h used a Fabry-Perot interferometer to study the velocity distribution 
of the low densi ty gas in the cavity and its relation to the cluster. A subsequent analysis 
bv lSmithI (119681) " suggested a n initial streaming velo city of 13 km s ^ at the ionization 
front. Later, in the study of [Fountain et all (Il97^ . 700 positions in the nebula were 



analyzed with multi-slit echelle spectroscopy in Hq, and these data revealed large sys- 
tematic velocities of up to 20 km-s~^ at the cavity, a sign of r apid evacuatio n of material 
due to powerful stellar winds. A similar result was found by 'Smith' ('1973'). 

The cluster itself has been the subject of many detailed studies over the years. The 
distance to this young cluster (and therefore to the entire complex) has been estimated 
many times with slightly different results. Table 1 is a compilation of these values, from 
which the most commonly used is 1600 to 1700 pc. 



O ptical Stud ies Some of the first visual photometric studies of NGC 2244 were made 

by iKirilloval ( I958h . who constructed the first HR diagram of the cluster and assigned 
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to it a distance modulus of 12.1 mag, and by I Johnsoiil (1962), who estimated the mean 
color excess in the clu ster to he E(B - V) = 0.46 for R = Ay/E jB - V) = 3.0. 
This was confirmed by iTumer J 1976h and lat er Qgura & Ishidal (Il98lh . who suggested 
a value of R = 3.2 it 0.15. lOgura & Ishidal (Il98ll) also propo sed an age of 4±l Myr 
and a star formation efficiency of 22% for the cluster. Later, iMarschall et all (1982) 
completed a proper motion study of 287 stars in the N GC 2244 area. They co nfirmed 
membership for 1 13 objects, 52 of them from the list of lOgura & Ishidal (Il98lh . 

A study t hat combined photometry as well as spectroscopy was completed by 
iPerez et ahl (11987.) . They found that some members of NGC 2244 presented anoma- 
lous values of R, possibly suggesting the coexistence of main sequence stars with very 
young objects-lik ely T Tauri stars. This was confirmed in the uvbyfi photometry study 
( Perez etaTIl 19891) . in which 4 mern bers presented evidence of being true pre-main se- 
quence (PMS) objects. Pe rez et all ([l989) also confirmed the age of NGC 2244 to be 
below 4 Myr but spread towards younger values, thus confirming a model of continuous 
formation. 
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Figure 3. An Ha vs. V-I diagram for NGC 2244. The solid line represents a 
ZAMS relation while the dashed line is a selection limit. Filled triangles are PMS 
stars while open triangles are PMS candidates. Bright members are marked with 
dark filled circles. X symbols are X-ray sources and dots are non-members. From 
[Park & Sung. C2002) . 



A study of great importance was performed by iPark & Sund (l2002l) . They ob- 
tained UBVI and Ha photometry for the cluster. They were able to determine member- 
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ship for a total of 30 cluster sources and to extend the hst of known PMS candidates to 
21. They subsequently identified members coincident with ROSAT point source cata- 
logs and compared their properties to the rest of the sample. S ix of the PMS candidates 
were confirmed as X-ray sources. In Figure [3l we show the iPark & Sund diagram of 
Ha emission vs. V-I color for NGC 2244. In this figure, PMS stars are clearly located 

above the main sequence. 

Using evolutionary models, iPark & SungI showed that most of the PMS stars and 
PMS candidates in their sample appear to have masses close to 1 Mq and an approxi- 
mate mean age of 0.4 to 0.9 Myr. Later, they estimated the main sequence turn-off age 
of the cluster to be 1.9 Myr, thus suggesting that the cluster has not stopped forming 
stars yet. 

Another important determination by Park & Sung is the Initial Mass Function 
(IMF) of NGC 2244. They found it has a flat (F = -0.7) IMF slope in the range 
0.5 < log 777 , < 2 .0, which compares well with the T = —0.8 slope o btained by 
Massey et all ( 1995h . By comparing directly to the IMF model of IScald ([l986) and 



to the observed mass function of NGC 2264, Park & Sung suggested that the shape of 
the stellar mass spectrum of NGC 2244 was dominated by high mass stars. However, 
they cautioned about the lac k of sensitiv i ty at t he intermediate and low mass ranges in 
their sample. More recently, IWang et all (l2008h . showed that NGC 2244 actually has a 
standard IMF shape regardless of the large number of massive stars. 



Spectros copic Studies The most complete spectroscopic study of NGC 2244 was 
done b y IVerschuerenI (119911). and it has bee n widely used in the literature. In par- 
ticular. Park & Sung used Verschuereiil (Il99ll) data to identify the spectral types of can- 
didate T Tauri stars in NGC 2244. Most of a vailab le spectral classifications for stars in 
NGC 2244 are com pile d bylO gura & Ishida (Il98lh . b ut additiona l data c an be found in 
iPerez et all (Il987h and iMassey et all (i 19951) (see also IChen et all |2004| and references 
therein). 

A low resolution, single slit investigation by iHensberge et all ( 19981) of 2 mem- 
bers and 3 field stars in the region of NGC 2244 yielded evidence that these were 
chemically peculiar, possibly magnetic stars. Hensberge et al. (2000) performed spec- 
troscopic analysis of the binary member V578 Mon, which resulted in an estimated 
distance slightly lower than other photometric estimates (see Table 1). The y also cal- 
culate d the age of the system to be 2.3±0.2 Myr. Another study by .Bagnulo et al. 
(|2004|) revealed that the source NGC2244-334 (a He-weak B5 V type star) has a very 
large longitudinal magnetic field (—9 kG), which is the second largest found in a non- 
degenerate star. They suggest that given the age of NGC 2244, this source might also be 
among the youngest Ap-Bp stars ever obser ved, sugg esting that these kind of sources 
can show magnetic fields from their birth. Li et al. ({2002) presented low resolution 
spectra for a sample of X-r ay counterparts from the ROSAT PSPC survey (see also 
Gregorio-Hetem et al.lll998h . They were able to confirm that five sources had strong 



Ha emission. Two of the stars were found to be Herbig Ae/Be and two others appear to 
be WTTS. These data indicate that X-rays are an efficient tracer of young populations. 



Near Infrared Recent surveys in the nea r-infra red have permitted investigation of the 
extension and structure of NGC 2244. Q (l2005h analyzed 2MASS data and suggested 
that NGC 2244 had a second component locate d approximately 6.6 pc west of the core 
center. Data from the FLAMINGOS Survey dRoman-Zuniga et all l2008h confirmed 
the existence of this second association, which is coincident with the area originally 
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Table 2. OB Members of the Young Cluster NGC 2244 



Star ID" RA*" DEC'' Sp. Type'^ K Other' 

noDD — — 



HD 46223 


06 


32 


09 


32 


-H04 49 


24 


6 


04V((f)) 


7 


32 


6 


68 


OI8 1-203 


HD 46150 


06 


31 


55 


52 


-H04 56 


34 


3 


05V((f)) 


6 


75 


6 


44 


OI81-122 


HD 46485 


06 


33 


50 


95 


-H04 31 


31 


6 


07V 


8 


20 


7 


45 


OI81-387 


HD 46056 


06 


31 


20 


87 


-H04 50 


03 


9 


08V((f)) 


8 


16 


7 


82 


OI81-84 


HD 46149 


06 


31 


52 


54 


-^05 01 


59 


1 


08.5V((f)) 


7 


59 


7 


25 


PS-160 


HD 258691 


06 


30 


33 


31 


-H04 41 


27 


6 


09V((f)) 


9 


70 


7 


93 


OI8 1-376 


HD 46202 


06 


32 


10 


48 


-^04 57 


59 


7 


09V((f)) 


8 


20 


7 


72 


PS-305 


HD 259238 


06 


32 


18 


22 


-^05 03 


21 


7 


BOV 


11 


10 


10 


28 


PS-348 


HD 46106 


06 


31 


38 


40 


-^05 01 


36 


3 


B0.2V 


7 


95 


7 


62 


PS-75;OI81-115 


MJD95 


06 


31 


37 


08 


-H04 45 


53 


7 


B0.5V 


15 


15 


12 


20 


— 


HD 259135 


06 


32 


00 


61 


-H04 52 


41 





B0.5V 


8 


54 


8 


12 


PS-226; OI8 1-200 


GSC 00154-00234 


06 


33 


37 


49 


-H04 48 


47 





B0.5V 


11 


89 


8 


64 


— 


HD 259012 


06 


31 


33 


46 


-H04 50 


39 


7 


BIV 


9 


35 


8 


79 


PS-43; OI81-80 


HD 259105 


06 


31 


52 


00 


-^04 55 


57 


3 


BIV 


9 


42 


8 


95 


PS-155; OI81-128 


GSC 00154-02337 


06 


32 


06 


13 


-H04 52 


15 


3 


BlIII 


9 


73 


9 


39 


PS-269; OI8 1-201 




Uo 


JJ 


J4 


A 1 
41 


_i_n/i "^Q 
+U4 jy 


A A 
44 


O 


xS 1 V 


/ 


OJ 


D 


o / 


r^TC 1 ^CQ 

*Ji?5t-j?5V 


BD-l-05 128 IB 


06 


31 


58 


93 


-H04 55 


39 


9 


B1.5V 


10 


38 


9 


74 


PS-214; OI81-193 


HD 259172 


06 


32 


02 


59 


-^05 05 


08 


6 


B2V 


10 


71 


10 


09 


PS-240; OI8 1-167 


GSC 00154-02247 


06 


33 


06 


56 


-^05 06 


03 


4 


B2 


12 


85 


11 


20 


OI8 1-345 


GSC 00154-02504 


06 


31 


31 


47 


-^04 50 


59 


6 


B2.5V 


10 


64 


10 


24 


PS-34; 0181-79 


GSC 00154-02141 


06 


31 


47 


89 


-^04 54 


18 


1 


B2.5V 


11 


66 


10 


87 


PS-123; OI81-130 


GSC 00154-02187 


06 


31 


58 


91 


-H04 56 


16 


2 


B2.5Vn 


11 


26 


10 


54 


PS-213; OI81-190 


GSC 00154-01007 


06 


32 


09 


84 


-H05 02 


13 


4 


B2.5V 


11 


22 


10 


43 


PS-300; OI8 1-172 


GSC 00154-01016 


06 


32 


24 


24 


-H04 47 


03 


7 


B2.5V 


11 


41 


10 


58 


PS-398; OI81-274 


GSC 00154-01247 


06 


33 


50 


56 


-^05 01 


37 


6 


B2.5V 


11 


12 


9 


99 


OI8 1-392 


GSC 00154-01753 


06 


32 


15 


49 


-H04 55 


20 


4 


B3 


12 


01 


10 


96 


PS-327; OI81-194 


MJD95 


06 


32 


22 


49 


-^04 55 


34 


2 


B3V 


15 


39 


13 


48 


PS-383 


HD 259268 


06 


32 


23 


04 


-H05 02 


45 


7 


B3 


11 


09 


10 


33 


PS-390 


HD 259300 


06 


32 


29 


39 


-^04 56 


56 


1 


B3Vp 


10 


79 


9 


34 


PS-441 


MJD95 


06 


33 


10 


16 


-H04 59 


49 


9 


B3V 


14 


98 


12 


39 




GSC 00154-02164 


06 


32 


51 


79 


-H04 47 


16 


2 


B5V 


12 


88 


11 


53 


OI8 1-334 



" Star identification. HD catalog is used when available, followed in priority by GSC an d BD 
catalogs. Stars listed as MJD95 are identified by position in the list of IMassev et al.l il995n . 
Coordinates from the 2MASS Point Source Catalog 



" Spectral types as listed bv lWang et al.l d2008h . except GSC 00154-02164 

Visual (V) photometry from 'Ogura & Ishida' ('1981"), Massey et al (1995) and iPark & Sund 
([2002). Near-infrared (K) photometry from the 2MASS Point Sour ce Catalog. 
" Correspondence with the lists of lOgura & Ishidal ( 119811) [OI] and lPark & Surig d2002h [PS]. 



labeled as NGC 2237. In Figure |4] we show the locations of infrared excess (IRX) 
sources with local surface densities above the field levels. The extinction in the region 
is rather low i{Av) ~ 3.0) compared to the molecular cloud areas but the distribution 
of excess sources is well defined and trace the location and extension of NGC 2244 and 
the satellite cluster. 

Using a cross-referenced catalog of NGC 2244 (Roman-Zuniga & Lada, in prepa- 
ration) which includes data from the 2MASS and FLAMINGOS near-infrared sur- 
veys and the survey of iPark & SungI (120021) . we were able to construct two optical- 
infrared photometric diagrams: Figure [5] is an optical-infrared color-magnitude dia- 
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Figure 4. Infrared Excess sources with significant local surface densities in the 
central part of the Rosette Nebula, from the FLAMINGOS survey. The clusters 
NGC 2244 and a possible companion cluster at the west region of the nebula, 
NGC 2237, are traced by this population. The back ground levels indicate 60 ii m 
emission from IRAS, in steps of 5.0 Jy. Adapted from lRoman-Zuniga et al.l ( l2008l) . 



gram in which probable members (P=0.6 or more) from the study of iMarschall et al. 



1982h are shown with dark symbols. The 1,2,3 and 5 Myr isqchron es from the pre- 



main sequence evolutionary models of 'P'A ntona & Mazzitellil ^199 1^ are also shown. 
The 1-3 Myr isochrones are in good agreement with our data (at least for tho se stars 
which are less affected by reddening). Notice that the lD'Antona & Mazzitellil models 
are calculated only for stars with 3 M© or less. 

Figure |6] is an infrared color-color diagram of NGC 2244. The scarcity of sources 
along the reddening band shows that the extinction towards NGC 2244 is rather low, 
as the parental material has been almost completely evacuated. A few of the compiled 
stars have large H-K colors and are located to the right of the reddening band with 



similar slope and direction as the classical T Tauri locus (iMever et al.lll997D, but with 



slightly bluer J-H colors likely caused by nebular gas extinction ( Heraudeau et al.ll 19961 : 
llshii et alll2002h . Another interesting feature is that some of the cluster members (solid 
symbols) form a locus that protr udes to the right of the reddening band, coincide nt with 
the observed locus for Be stars dHowells et al.ll200"ll : ISubramaniam et al I2005h ,which 
suggests the presence of circumstellar material around some of the massive members 
of the association. 

In Table |2] we present a list of the O and B members of NGC 2244, as listed in 
the X-ray studies by Townsley et al. (2003) and Wang et al. (2007). The list includes 
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NCC 2244 




1-1 [r,.,y] 



Figure 5. An optical-infrared color-magnitude diagram for NGC 2244. The solid 
diamond symbols have membership probabilities higher than 60%. The broken lines 
are model PMS isochrones at 1, 2, 3 and 5 Myr from D'Antona and Mazzitelli 
(1997). The thick sohd line is the ZAMS. 



identifications, positions, spectral type s, B and K photortietry an d a cross identification 
with the lists of lOgura & Ishidal dlQSlh and lPark & Sungl (1200211 . 



Mid-Infrared Observations NGC 2244 was observed with the Infrared Array Cam- 
era (IRAC) and the Multiband Imaging Ph otometer (M I PS) o n board the Spitzer Space 
Telescope in 2004 and 2005, respectively. iBalog et aP (l2007l ) obtained from these ob- 
servations a catalog with more than a thousand high quality sources detected simultane- 
ously at the four IRAC bands (3.6, 4.5, 5.8 and 8.0 /zm), plus a total of 279 MIPS coun- 
terparts at 24 /zm, succeeding in completing a census of NGC 2244 in the mid-infrared 
down to a mass limit of --^0.8 Mq. The main purpose of the study of iBalog et al.i was 
the investigation of circumstellar disk survival in the OB environment. 

A total of 337 and 25 sources were successfully identified as Class II and Class I 
sources, respectively (with 213 and 20 sources identified likewise in the IRAC-MIPS 
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Figure 6. A near-infrared color-color diagram of NGC 2244. The solid diamond 
symbols have membership probabilities higher than 60%. The solid black lines indi- 
cate the ZAMS for dwarfs and giants, while the dotted lines indicate the reddening 
band from the law of Cohen (1981). 



combined catalog). This implies that the ratio of Class I to Class II sources in NGC 
2244 is less than 7%, much lower than t he ratios obs erved in star forming regions 
with smaller proportions of massive stars. iBalog et all found that Class I sources are 
preferentially located at the outer regions of the cluster, with none located near the 
central parsec, opposite from the Class II sources which appear to concentrate around 
the OB center. These two aspects suggest that within the cluster star formation might 
have stopped, but it continues outside the HII region (see Figure |71l. 

Later, Balog et al. divided the detections into stars with and without infrared 
excess to give an estimate of the disk ratio in the cluster. They found that the disk ratio 
within 0.5 pc from the center of the cluster is 27%, then it increases to about 45% within 
0.5 and 2.5 pc. Th e average disk fraction is thus above 40%, in close agreement with 
Haisch et all d 200ll) which suggest that the disk fraction in a cluster decreases from 80 
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Figure 7. Distribution of candidate Class 1 sources (diamonds). Class 11 sources 
(circles), Class 1-11 sources (squares) and O stars (solid circles) in NGC 2244, over- 
laid on a 24 /im image obtained with MIPS onboard the Spitzer Space Telescope. 
Adapted from .Balog et al.. (.2007.) . 

to 50% within the first 2-3 Myr of life of a clusteiQ. The disk fraction in NGC 2244 
calculated by Balog et al. is not very different (A ?s 10%) from those found in clusters 
of similar age without O stars, plus they found a steep decrease of the photoevaporation 
rate as a function of distance from the cluster center, allowing them to conclude that the 
effect of high mass stars on disks is only significant at very close distances to the hot 
stars. 

X-Ray Studies NGC 2244 is an important target for X-ray studies due to the interest 
in investigating the nature of massive stars as source s of hig h energy photons. Ob- 
servations made with the Einstein space telescope by iLeahyl (fl985) revealed that the 
OB stars in NGC 2244 were embedded in a bubble of hot, low-density gas with dif- 
fuse X-ray emission at a temperature of 1 keV, possibly originated by hot shocks in 



For comparison, the disk fraction observed in the near-infrared is 12% (see Section 4.1), but such dis- 
crepancy is due to the fact that IRAC is more successful at detecting circumstellar material because the 
cold and hot dust that comprises circumstellar disks emits preferentially within 3 and 10 fim 
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the stellar winds of O stars. The ROSAT Consortium observations (I2000h yielded 34 
X-ray sources in NGC 2244, with typical X-ray luminosities of IQ^O _ iq32 ^^^^ 
Six of these X-ray sources a re PMS candidates as reported by iPark & Sun 3 (l2002h . 
iBerghofer & ChristianI (l2002h also studied NGC 2244 ROSAT sources and found that 
objects with the faintest X-ray emission have very high X-ray to optical luminosity 
ratios. They noted that the number of X-ray emitters associated with Ha emission in 
NGC 2244 is remarkable. 
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Figure 8. A map of the Rosette Complex. The background image is a plate of 
the IRAS survey at 25 um. The con t ours r epresent ^^CO levels (as Tp^ak) from the 
survey of ^ Hever. Williams. & Brunl (l2006h . Labels indicate the main regions of the 
molecular complex, originally identified by Blitz & Thaddeus (1980). 



A more recent study by IChen et all ( 20041) was dedicated also to analyze the co- 



incidence of ROSAT with stars of NGC 2244. A careful re-samphng of the ROSAT 
data combined with new optical imaging and spectroscopic data yield 47 X-ray peaks 
with stellar counterparts in the Nebula region. Although it was estimated that about 
half of the counterparts were late type foreground stars, they could confirmed X-ray 
emission for the rest. Specifically, X-ray emission was confirmed for the 2 earliest type 
stars in the cluster (HD 46223 and HD 46150 with types 05 and 04 respectively), for 
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a very young binary (J063148.30+045820.5), for the s ource HP 259210 \ yhich coin- 
cides with the location of high speed knots identified by lMeabum & WalshI ([1986'), and 
fo r the source J063309 .61+044624.5 located near a massive nebular pillar identified 
by ISchneps et alj (119801) (see also Section 5.2). Taken together, all these investigations 
give strength to the hypothesis that X-ray emitters in NGC 2244 are either very massive 
or very young in the case of late types. 

Probably the most sensitive X-ray observations of NGC 2244 to date were ob- 
tained with the Advanced CCD Imaging Spectromete r (ACIS) on board the Chandra 
space Observat ory. The point source a nalysis by Wang et al.l (120081) is complementary 
to the study of iTownsley et all (l2003h which was dedicated to inves tigate the di ffuse 
emission in the HII region (see Section 3.2). Within the ACIS mosaic. 'Wang et all were 
able to extract, with high reliability, over 800 X-ray point sources within 20 arcmin 
from (a, (5) = (6''31™59^ +4'='55'36", J2000). The X-ray sources have luminosities 
between 10^^ '' and lO^^.o 

ergs-s ^ in the hard band (2-8 keV). All sources with spectral 
types earlier than B 1 were detected. 

A total of 7 12 Chandra sources in the lWang et all survey were found to be counter- 
parts of near-infrai^ed sources detected in the FLAMINGOS (see Section 4) and 2MASS 
surveys, down to a mass limit close to 0.1 Mq. Among this X-ray selected sample, the 
fraction of stars with mass > 0.5 with significant K-band excess (likely indicative 
of a circumstellar disk) is about 6%, and three of the sources were positively identi- 
fied as potential Class I sources. The unobscured population in NGC 2244 is about 
1.2 times larger than that of the ONC (Feigelson et al. 2005) and both the XLF and the 
KLF (constructed from the uniformly sampled 2MASS catalog) appears to suggest a 
normal Salpeter IMF for stars with mass > 0.5 Mq. By comparing to the ONC dis- 
tribution, they estimate that the total population of NGC 2244 could be around 2000 
members. The radial source density profile of X-ray sources appears to have a relaxed 
structure around the center of the cluster and suggests a low probability of significant 
mass segregation. 



3. The Rosette Molecular Cloud and Embedded Populations 

3.1. Structure of the Rosette Molecular Cloud: CO studies 

First attempts to detect CO emission associated with the Rosette nebula were unsuc- 
cessful as they pointed at the nebula region, which is mostly composed of neutral and 
ionized gas. The observations reported by Blitz & Thaddeust (il98Q) . which targeted the 
southeast adjacent region of the nebula, were the first successful detections of molec- 
ular gas in the Rosette. Their NRAO survey mapped over 80% of the ^^CO emission 
in the area of the cloud with a 1' beam size, and yielded information about its large 
scale distribution. They estimated the angular extent of the cloud to be 3.5° (98 pc at a 
distance of 1600 pc) and labeled the most prominent sub-structures (see Figure[8]l. 

Particularly important regions are the Monoceros Ridge (region A 1-2) which is 
literally a region of gas compression at the cloud-nebula interface; the central core 
(Al-1), which hosts the most massive clumps in the cloud and is the strongest region 
of star formation; the cores D and G, which are separated from the main body of the 
cloud but have ongoing star cluster formation; the IRS core, which hosts the massive 
proto-binary AFGL-961 (see Sections 4 and 5.1); the back core B, which is more loose 
in structure than the regions near the nebula; and the arm or E core, which despite 
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its brightness contains no significant star formation (no IRAS sources, or near-infrared 

clusters have been found in this core so far ). 

In a subsequent study, Blitz & StarkI (1 19861) mapped the ^^CO and ^^CO emis- 
sion with improved sens itivity at the AT&T Bell Labs, unc overing the high degree of 
dumpiness of the cloud. IWiLIiams. de Geus. & Bhtzl {1994) made use of the data from 
iBlitz & Stark and listed a total of 95 clumps. The clumps with evidence of star for- 
mation had larger peak temperatures, larger densities and als o were more gravitation- 
ally b ound compared to clumps from the Maddalena Complex ([M addale na & Thaddeus 



19851) . a cloud with very low star formation. Later, IWilliams et al. (19951) showed that 



about half of the clumps in the Rosette Molecular Cloud were gravitationally bound and 
the rest were supported by pressure from the inter-clump medium, which was shown to 
be mostly atomic and about 40 ti mes less dense. In Fig ure [9| we show the locations and 
relative sizes of the clumps from lWilliams et all (1 19951) . 
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Figur e 9. Lo cations and relative sizes of clumps in the Rosette Molecular Cloud 
from I Williams et al. (1995). The size of the symbol is proportional t o the mass of 
the clu mp. The contours represent CO emission from the survey of iBlitz & StarkI 
and the solid line delimits their coverage in "CO. 



From the clump central velocities IWilliams et al.l (|1995h found that the cloud has 
a well defined velocity gradient of about 0.08 km-s~^-pc~^. Also, the negative corre- 
lation between clump mass and clump to clump velocity dispersion sugg ested that the 
system is still far from equipartition even though it is dynamically evolved. IWilliams et al 
(1 1995b also found the star forming activity to be more intense in the ridge and core ar- 
eas, near the interface cloud-nebula: clumps located near the nebula presented larger 
excitation temperatures, average densities and star forming efficiencies, all clues of 
cloud evolution. Other properties of the clumps (mass, sizes and line widths) did not 
show any significant variations along the cloud. 
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Another CO study was done by lSchneider et all (Il998al) . The observations focused 
on the central part of the cloud, detailing the structure of the midplane star forming 
cores. They paid special attention to the IRS core, where the source AFGL-961 is 
located (see also Section 5.1) and pointed out the ample blue wing emission due to the 
powerful outflow from this object. In a complementary study, ISchneider et al.l ( 1998bh 
examined the CII emission (158 ^m) at the ridge, the central core and the IRS core. 
They found weak but significant C+ emission deep into the molecular cloud cores and 
suggested that the distribution agrees well with a clumpy molecular cloud exposed to a 
low level UV radiation field. The penetration of UV photons in the cloud is apparently 
facilitated by a high density contrast clump-interclump medium. 

The clump mass spectrum in the Rosette Molecular Cloud has the form dN/ dM oc 
M^^ where x ^ 1.6, depending on the range, bin size a nd beam res olution used. 
The exponent in this power law is similar to other clouds (lBlitz|[l993 ). but what is 
more important, it is much shallower than the one corresponding to the stellar IMF 
(x = 2.35). This shows that although there are more small clumps by number, most 
of the mass is contained in only a few big clumps, while for stars both numbers and 
total mass are dominated by the lowest mass bins. Interestingly enough, the power 
law index is in fact similar to that corresponding to the mass distribution function for 
embedded clusters, which is suggestive of a uniform star formation efficiency for most 



star forming cores ( Lada & Ladal 



The more recent survey of 



2003h. 



Heyer. Williams. & Bruni (l2006h . obtained with the 



wide field array SEQUOIA at the FCRAO I4m telescope achieved resolutions of 45" 
at 115 GHz and 47" at 110 GHz. The maps revealed "textural variations" in the ^^CO 
emission across the complex, with a brighter emission component within the nebula 
projected radius (approx. 30 pc from the center as de fined by Ce lmk, 1985) and weaker , 
more extended emission outside this ionization edge. iHeyer. Williams. & Brunt (I2006h 
suggest that the weaker emission is probably due to subthermally excited material with 
lower densities. They also calculated the total molecular mass of the cloud to be 1.6 x 
10^ M0 from ^^CO, and found a LTE mass of 1.16 x 10^ Mp, from ^^CO. Moreov er, 
they were able to apply a Principal Component Analysis (Hever & Schloerbl ll997 ) to 
determine the turbulent flows and the turbulence scale in the Rosette Molecular Cloud. 
The analysis revealed more significant variations in the velocity structure of the cloud at 
the regions located within the ionization than in the more diffuse, external component. 
This fact confirms the interaction of the cloud and the HII region. They suggested, 
however, that the interaction is still very localized, and have not affected the global 
dynamics of the cloud yet. 



3.2. Interaction with the Rosette Nebula 



Celnikl(ll986h focused on comparing his Ha map and radio continuum obser vations of 



the neb ula (see Section 2.1) with the CO map of the molecular cloud from Blit z & Thaddeus 



(11980). Celnik constructed a complex model of the distribution of the main CO cores 
(see Fig 121) in the context of the HII region and estimated the rotation center of the cloud 
at (a, 6) = (6''32'"39^ +4'^19'43", J2000). Finally, he re-calculated the mass of the 
entire complex by adding the total mass of ionized atoms, stars, dust and molecular gas, 
resulting in 3.3 x 10^ Mq. 

Another important study was done by lCox et"aL who made use of the avail- 



able infrared emission data from the IRAS satellite (12, 25, 60 and 100 /um). ICox et al" 



were able to study in great detail the distribution of dust and compared this to the dis- 
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tributions of ionized and molecular gas. Additionally, they were able to estimate a total 
infrared luminosity of roughly 1.1 x 10^ Lq, or about 50% of the available luminosity 
from the cluster NGC 2244. Warm dust (usually present near an O B association) typ- 
ically emits strongly at the four IRAS bands. However, ICox et all showed that in the 
Rosette, while the 60 and 100 ^m emission are quite strong at regions of ionized and 
neutral gas (nebula), the 12 fim emission is preferentially located beyond the limits of 
the ionization front (molecular cloud), suggesting a heavy rate of destruction of dust 
grains from UV radiation from the cluster. Surprisingly, the 25 ^m emission was found 
to be significant in some parts of the ionized nebula, possibly due to the existence of a 
second ty pe of dust particle that is more resistant to UV photons. This was also sug- 
gested by IShipman & Clark* (1994), who found that the maximum temperature in the 
shallow temperature gradient found in the nebula (see Section 2.1) seems too low to 
sublimate ice mantles in grains and too low for grains to emit significantly at 12 or 
25 Mm — a proble m possibly solved with a second type of grain in the region. Later, 
Shipman & Careyl (11996) suggested that line emission could be contributing strongly 
to the IR emission of the nebula, and suggested once more that the presence of a "hot 
dust" component is necessary to model this emission, especially at 25 fim. 

Figure [To] shows the superposition of the 12 fim emi ssion from the IRA S survey 
and a 1400 MHz radio continuum emission map from Hol daway et al. (Il999h . The in- 
frared contours indicate that the warm dust emission defines a shell that encloses the 
ionization front, showing the effect of heavy dust destruction by the nebula. The super- 
position of these maps defines very clearly the region where the HII bubble impacts the 
molecular cloud. 



Patel et all (119931) observed the SE quadrant of the nebula, which coincides with 



the regions located south and ea st of the Monoceros Ridg e, where molecular gl obules 
were identified by iBlockl (1 19901) . iBlock et all (Il992h and ISugitani et all (Il99lh . The 
CO and ^^CO (J = 1 ^ 0) data of iPatel et all confirmed that this part of the cloud 
contained nine well defined cometary globules, which appear to be blue-shifted by ap- 
proximately 6 km-s^^ with respect to the Rosette Molecular Cloud. The globules have 
sizes of 1-3 pc and masses of 50 to 300 Mq and four of them are coincident with IRAS 
sources and local maxima in the ^^CO emission. Patel et al.i suggest that these globules 
cannot be formed as Rayleigh-Taylor instabilities as in the case of globules (elephant 
trunks) in the NW quadrant of the Rosette Nebula (see Section 5), but rather by a mech- 
anism of radiation driven implosion, caused by penetration of UV radiation from the 
nebula cluster which favors D-type ionization front s inside mas s ive cl umps, causing 
the formation of the globules. A follow up study by 'White et ajj ( 1997h also included 
CO, ^^CO, and C^^O observations of a similar region, revealing that the Globule 1 from 



Patel et aLl is indeed the largest cometary globule in the area; this object has a prominent 
head-tail structure with a projected length of almost 1.3 pc which points away from the 
center of the Nebula. They also confirmed the agreement with a RDI model and were 
able to calculate an age of 0.4 Myr for the globule, suggesting that it has almost reached 
the stage of maximum compression and is close to a quasi-static cometary phase. 



High Energy Studi es When the Rosette Mole cular Cloud was confirmed as a region 
of star formation, Gregorio -Hetem et all (Il998h used ROSAT data again, this time to 
map the MonR2 cluster and the Rosette Molecular Cloud areas in order to confirm 
a correlation between star forming cores and clusters of X-rays sources. They found 
strong X-ray emission in NGC 2244, at the ridge of the cloud (A 1-2 in fig[8]l, and at the 
cloud core area (Al-1), but the resolution was poor and individual sources could not 
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Figure 10. IRAS 12 /im e mission map supe rimposed on a 1400 MHz radio con- 
tinuum emission map. Credit: iHoldaway et alJ (1999) and NRAO. 



be resolved. They suggested that molecular cores known to have active star formation 
but failing to show significant X-ray emission, could be predominantly forming low- 
mass stars. They also suggested that detectable X-r ay counterpart s are in most case s 
Herbig Ae/Be or T Tauri stars, as found in NGC 2244 (iLi et al.ll2002l : IWang et all2008l). 

^The more recent observations of the Rosette Complex done with Chandra (Towns ley et al. 

I2OO3I) have resolutions of only a few arcseconds, thus allowing for the detection of X- 
ray counterparts for 75% of the OB members of NGC 2244. One of the most interest- 
ing results of this study was the confirmation of a second, soft diffuse emission which 
probably originates from the O star winds and is later brought to thermalization by 
wind-wind interactions or by the shock with the surroundings, in this case the molecu- 
lar cloud (see Figure [TTI). This X-ray plasma surrounds the OB association and fills the 
nebula cavity completely. 
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Figure 11. A 0.5-2 keV Chandra image of the Rosette Complex. The emission 
has been smoothed to highlight the soft diffuse emissio n that originates in the neb- 
ula and propagates into the molecular cloud. Credit: iTownsley et al.l (2003) and 
NASA/Chandra X-Ray Observatory. 



The Rosette Nebula may not only be interacting with its companion m olecular 
cloud. The location of the Rosette Nebula near the edge of the Monoceros Loop dPavies 
— also known as the Monoceros Supernova Remnant — at the edge of the NE 
quadrant, has motivated a number of studies aimed at investigating the high energy pho - 
ton emission in that region. Deep Ha-i-[NII] photographic plates by lDavies et al. ( 1978 ) 
suggested a correlation between a filamentary structure observable in Ha emission and 
a Rosette nebula feature observable in 240 MHz radio waves. This featu re was proposed 
as evidence of loop-nebular interaction and confirmed by decameter (Odegard 198^ 
and diffuse X-ray emission (Leahy et al. 1986) observations. Later, high energy (100 
MeV) 7-ray images from EGRET (Jaffe et al. ,1997i) . revealed a feature partly coinci- 
dent with the filaments and apparently significant (la) over expected diffuse emission. 
If real, this 7-ray feature would possibly be due to the interaction of charged particles 
with the dense ambient medium at the shock region. The HEGRA system of atmo- 
spheric Cerenkov telescopes at La Palma Observatory was used to calculate the cos- 
mic ray emission from the loop-nebula interaction region, and found one possible TeV 



emission source (3EG J0634-I-0521) with low statistical significance (lAharonian et al. 
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20041) . More recently, the HESS telescope array in Namibia was used again to map the 



SNR-Nebula interaction region and found one gamm a ray source, HE SS J0632+057, 
located close to the rim of the Monoceros Loop CAharonian et al.ll2007h . 



3.3. The Rosette Molecular Cloud: Embedded Populations 

Near-infrared Observations The coincidence of massive clumps and luminous IRAS 
sources pointed out by the study of IWilliams et all ( 1995h strongly suggested that star 



formation had already taken place across the molecular cloud. However, the poor spa- 
tial resolution of the IRAS point source survey did not allow the resolution of Individ 



ual rt iembers of an embedded population. Early near-infrared studies (e.g. iPerez et al. 
Il987h did not cover the molecular cloud areas, and optical photometric studies were 
incapable of detecting obscured populations. 



An exploratory near-infrared survey (JHK) by iPhelps & Ladal (1997) that made 



use of the imager SQIID finally confirmed the exi stence of embedded clusters in some 



of the most m assive clumps fr om the list of Williams et aLl that were associated with an 



IRAS source. iPhelps & Ladal were able to identify by visual inspection of their images 
seven deeply embedded clusters with bright nebulosities. Other massive clumps did not 
contained a cluster, which suggests that high mass is a necessary but not a sufficient 
condition for cluster formation. 

Complete spatial coverage of the Rosette Complex in the near-infrared was first 
accomplished with the release of the All-Sky 2MASS survey catalogs. The 2MASS 
survey was a major gain in data uniformity but unfortunately not in sensitivity. Due to 
the distance to the Rosette (1.6±0.2 kpc), the 2MASS comp leteness hmi t (K=14 .3 mag) 



is not deep enough to study the low mass end of the IMF. iLi & SmithI (120051) studied 



the distribution of 2MASS sources with large H-K colors and they were able to suggest 
three general areas of cluster formation. The extensions of these star fo rmation regions 
were defined visually, and could not be compared with the clusters of iPhelps & Ladal 



([1997) because at that point the extents and membership of individual clusters were not 

determined. 

The survey of iRoman-Zuniga et all (1200 8*) aimed for deeper observations using 



FLAMINGOS on the KPNO-2.1m. The project was done as part of a FLAMIN- 
GOS/NOAO survey of Giant Molecular Clouds Q. Observations of the Rosette covered 
the regions surveyed by 'Phelps & Lada' and extended to the Back Core and the neb- 
ula. Identification of young embedded populations was done by selecting objects with 
near-infrared excess characteristic of circumstellar emission. Also, a Nearest Neigh- 
bors analysis was applied to this color selected sample to facilitate the detection of 
embedded clusters as regions with significant local surface densities. 

The technique allowed for a cl ear confirmation of the location and extensions of 



the seven bright clusters reported in Phelps & Ladal labeled as PLOl to PL07, but also 



indicated the existence of another two clusters associated with massive clumps in the 
molecular cloud, labeled as REFL08 and REFL09. The sources with near-infrared ex- 
cess also traced well the extension of NGC 2244 and a possible new cluster at the west- 
em edge of the HII bubble, located at the region historically identified as NGC 2237, 
(see also Section 2.2). Finally, a small cluster was detected at the northeast end of the 



^see |http://flainingos.astro.ufl.edu/sfsurvey/sfsurvey.html| 
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nebula and was labeled REFLIO. The location of the clusters and their distribution of 
near-infrared excess sources is shown in Figure [121 
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Figure 12. Identification of young clusters in the Rosette Complex. Plus symbols 
indicate the location of NIRX objects brighter than K< 15.75 and local surface den- 
sities higher than the background average. Black dots are sources with near-infrared 
excess colors and local surface densities below the background average. The contour 
levels indica te a baseline level (0.8 K km s^^) of ^"^CO emission from the survey of 
iHeyer. '\^lhams. & Brunt ( 2006) used to define the projected surface of the cloud. 
The thin solid line marks the coverage of the FLAMINGOS survey. 

The sizes (equivalent radii) of clusters in the Rosette are relatively large, appear 
to be anti-correlated with their mean average extinctions, and with the near-infrared 
excess fractions. If larger near-infrared excess fractions indicate slightly younger ages, 
then the anti-correlation might be indicative of rapid cluster evolution: clusters could 
form as compact units then expanded in a timescale shorter than the T Tauri phase. 
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From the direct counting of near-infrared excess stars, it was possible to determine 
the fraction of the young stellar population that is forming in clusters at the present 
day. The 9 largest clusters embedded in the Rosette Molecular Cloud occupy approx- 
imately 242 C)^ or approximately 9% of the main molecular cloud areas covered by 
the FLAMINGOS survey. After correcting for background stars the embedded clus- 
ters account for 86±5% of the embedded population in the cloud. These fractions are 
comparable to what was found in o ther nearby molecular clouds like Orion or Perseus 
(lLadaetalJll99ll : ICarpeDteill2000h . where clusters comprise 50 to 96% of the young 
stars that can be observed. 



Table 3. Young Clusters in the Rosette Complex 



ID" 




dec'' 


Clump" Rclump'^ 


No. IRX" 


No. Class I-'^ 






(J2000) 




[pc] 


K< 15.75 


M<0.4 Mq 


[pc] 


PLOl 


06 31 49.32 


+04 19 34.7 


11 


2.21 


29±5 


9 


1.16 


PL02 


06 33 16.78 


+04 35 32.1 


18 


1.33 


32±6 


12 


1.46 


PL03 


06 33 33.07 


+04 00 11.2 


7 


2.79 


80±9 


>6 


1.69 


PL04 


06 34 13.92 


+04 25 05.4 


1 


2.36 


89±9 


<107'' 


1.85 


PL05 


06 34 30.70 


+04 20 01.8 


19 


1.50 


57±8 


<107 


1.31 


PL06 


06 34 38.76 


+04 12 55.4 


2 


3.06 


13±4 


>10 


0.75 


PL07 


06 35 29.78 


+03 59 10.9 


3 


2.79 


22±5 


27 


0.88 


REFL08 


06 34 18.98 


+04 20 03.6 


17 


1.71 


49±7 


<107 


1.30 


REFL09 


06 35 07.73 


+03 41 34.7 


5 


3.60 


65±8 


12 


1.49 


PouC 


06:33:12.68 


+04:31:00.5 








12 


0.90 


PouD 


06:33:40.98 


+04:03:56.1 








9 


0.83 


REFLIO 


06 31 06.79 


+05 14 50.0 






15±4 




1.15 


NGC 2237 


06 31 58.51 


+04 54 35.7 






36±6 




1.91 


NGC 2244 


06 30 36.10 


+04 58 50.6 






62±8 


176 


2.30 



a ID key: PL— Phelps & Lada (1997); REEL— Roman-Zuniga et al. (2006a); Pou— Poulton et 
al. (2008) 

''Location of peak surface density of near-infrared excess sources 



'^Molecular clump host identification from Williams, Blitz & Stark (1995) 
'^Clump host equivalent radius 
"Limited to K< 15.75 mag. 

■'^Adapted from Table 2 and figure 9 of Poulton et al. (2008). Mass detection limit is 0.4 Mq 

^Cluster radii defined from IRX local surface density contours, except PouC and PouD, listed 
as equivalent radius. 

''Clusters PL04, PL05 and REEL08 are reported as a single cluster in Poulton et al. 



Mid-infrared Observations Recently, iPoulton et all (l2008h reported the first results of 
the Spitzer telescope survey of the Rosette Molecular Cloud. The area surveyed with 
IRAC and MIPS covers 1x1.5°, similar to the one studied with FLAMINGOS. The 
Spitzer observations allowed to identify, via the fit of SED profiles, a total of 75 1 young 
stellar objects with infrared excess down to a mass limit of 0.4 M0. Also, they used 
a method of nearest neighbors applied to all sourc es with infr ared excess to identify 
embedded clusters, confirming the seven clusters of iPhelps & L ada ( 1997) as the major 
clusters in the cloud and confirming the existence of clusters REFL08 and REFL09 
(labeled as clusters E and F, respectively). Their analysis also revealed two new small 
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clusters or associations, labeled clusters C and D (we used the labels PouC and PouD 
in this work). These two groups are located near clusters PL02 and PL03, respectively. 

The Spitzer survey of the Rosette confirmed that Clusters PLOl to PL05 have very 
dense envelopes of hot dust, evidenced by strong extended emission in all IRAC and 
MIPS bands. The Spitzer survey also revealed that sources within 15 pc from the center 
of NGC 2244 are definitely older than those deeply embedded in the cloud, with 2% vs 
7% of identified Class I sources, respectively. 

In Figures \T3\ [141 [13 and [16] we present false color images constructed with 
near-infrared images from FLAMINGOS and with mid-infrared images from IRAC 
for all the embedded clusters in the Rosette Molecular Cloud: PLOl-07, REFL08-09 
and PouC-D. The FLAMINGOS images, with higher resolution, define well the stellar 
content of the clusters, while the IRAC images show in addition the distribution of hot 
dust and gas around the embedded clusters. 

Table [3] contains information about the location of the 14 known clusters em- 
bedded in the Rosett e Molecular Cloud, along with the corresponding clump from 
WilUams et al. I (Il995h . Also indicated is the number of near-infrared excess stars de- 



tected in each cluster region and the cluster equivalent radii, which was calculated by 
estimating the extension of each cluster from local surface density contours of near- 
infrared excess sources above the backgr ound level. We also indicate the number of 
Class I sources reported by iPoulton et al.l 

Interaction of Embedded Populations with the Local Environment In order to give 
insight into the nature of the embedded populations of the Rosette Molecular Cloud, 
a high resolution mm- wave survey of clusters PLOl to REFL08 was performed at the 
IRAM 30m telescope (Roman-Zuniga 2006, Roman-Zuniga, Williams & Lada, 2008 
(in prep)). The cluster regions were scanned for emission in 6 molecules: ^^C0(2-l), 
^3CO(2-l), C^^Oil-l), CS(2-1), HCO+(1-0) and N2H+(l-0). All of the clusters show 
prominent clumps in most tracers, except cluster PL02 which has weak emission in 
C^^0(2-l) and N2H"'"(l-0). In Figure [TT] we show one example of the maps of dense 
gas emission maps for cluster PL06 (host of AFGL 961). 

Bulk physical properties of the clumps are calculated from these observations. 
Particularly, mass, sizes and velocity dispersions appear to be fairly similar for all of 
the cluster areas, which translates to no significant variation of the properties of star 
forming clumps across the cloud, as suggested by Williams et al. (1995) . However, the 
morphology of the gas emission is different from clump to clump, suggesting different 
stages of the clusters evolution. For example, while the gas emission around clusters 
PLOl, PL06, PL07 and REFL08 is compact and forms a well defined envelope around 
the clusters, clusters PL03, PL04 and PL05 appear to be partially emerging from the 
clumps, with the gas apparently being swept out from some areas of the cluster and 
recollected in others. 

The emergence of the massive clusters PL04 and PL05 might suggest a late stage 
of embedded cluster evolution. The emission of dense gas tracers for these clusters is 
confined to only two relatively small areas which coincide with the most obscured and 
dense parts of the respective cores. The rest of the stars in the clusters are still embedded 
in the larger and more quiescent ^^CO core but some of them are already visible in DSS 
optical plates. Clusters PL04 and PL05 are the largest and most extended embedded 
clusters in the Rosette Molecular Cloud (see Table[3l) but could also be the oldest. These 
clusters are located in the part of the cloud that is in direct interaction with the shock 
front from the Nebula, as described by iHeyer. Williams. & Brunt (l2006h . 



Figure 13. Embedded clusters PLOl, PL02 and PL03 (top, middle and bottom, 
respectively) in the Rosette Molecular Cloud. Left Panels: False 3 color images 
(J,H,K) constructed with images from the FLAMINGOS survey. Right Panels: False 
3 color images ([3.6], [4.5] and [8.0];um) constructed with images from the IRAC 
survey. AH images cover a 5'x5' area. Scale bars of 1.0 arcmin are also shown on 
top. 
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Figure 14. Same as figure[T3] for embedded clusters PL04, PL05 and PL06 (top, 
middle and bottom, respectively). Note: Artifacts in cluster PL06 IRAC image are 
due to saturation of IRAC detector by the bright source AFGL961. 



Figure 15. Same as figure[T3] for embedded clusters PL07, REFL08, and REFL09 
(top, middle and bottom, respectively). Note: IRAC picture for REFL09 was con- 
structed with [4.5], [5.0] and [8.0] /im, as it misses coverage in [3.6] /im. 
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Figure 16. Same as figure [13] for embedded clusters PouC and PouD (top, and 
bottom, respectively). Note: The area covered in these images is 3'x3'. 



The Sequential Formation Hypothesis In Figure [18] we show the distribution of the 

surface density of near-infrared excess sources (young stars) in the Rosette Complex as 

a func tion of the distance from the center of NGC 2244, as calculated by Roman-Zuiiiga et al, 
(120081) . Young clusters appear as peaks in this distribution, and it is possible to distin- 
guish four main groups of clusters. The most prominent group contains NGC 2244, 
NGC 2237 and REFLIO, which are the clusters located in the nebula. The second 
group contains clusters PLOl and PL02 and coincides with the cloud ridge (core Al- 
2). The third group contains clusters PL04, PL05, REFL08 and PL06, all located at 
the RMC central core, and cluster PL03 located in core D. The fourth and last group 
contains clusters PL07 and REFL09 which are located at the back core of the cloud. 
In the bottom panel of Figure [18] we show the average infrared excess fraction for each 
of these cluster groups. The average excess fraction appears to increase as a function 
of distance from NGC 2244, suggesting that the embedded clusters are progressively 
younger the further they are from the Nebula. 

The existence of this apparent age sequence indicates that star formation in the 
Rosette did indeed take p lace sequentially in time, but not in the way proposed by 
Ehnegreen & Ladal (Il977l) . because that model requires one episode of cluster forma- 
tion to trigger the next one. The relative age differences of clusters cannot be large 
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Cluster 6 
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Figure 17. Maps of dense molecular gas emission towards cluster PL06 in the 
Rosette Molecular Cloud, from Roman-Zuniga (2006). The first two panels show a 
false color near-infrared image of the cluster and a IRX surface density map. The 
molecular map panels have base and step values for the contour levels. The maps are 
integrated intensity except for ^^CO(2-l) which is mapped as peak temperature. 



enough to account for such a scenario. Instead, it appears that the overall age sequence 
of cluster formation may be primordial, possibly resulting from the formation and evo- 
lution of the molecular cloud itself. The HII region cannot be responsible for the se- 
quence of cluster ages, although it does appear to have a significant impact on the 
underlying sequence by either enhancing or inhibiting the star forming process: At the 
cloud ridge, the direct interaction with hot, ionized gas might contribute to the rapid 
evaporation of cluster envelopes in PLOl and PL02, possibly stopping the star forma- 
tion process at an early stag e; however, re cent Chandra and Spitzer telescope observa- 
tions (Wang et al. (in prep); Poulton et al.l . 2008) suggest that PL02 hosts two Class O/l 
sources, which is suggestive of cluster substructures younger than PL04 or PL05. At 
the central core, where the main interaction between the nebula and the molecular cloud 
takes place, the result is an enhancement of the star formation efficiency, evidenced by 
the fact that the clusters in this area (PL04, PL05, PL06, REFL08) account for almost 
50% of the total embedded population in the cloud. Finally, at the back core the influ- 
ence of the HII region is minimal and star formation in this part of the c loud may have 
occurred spontaneously dRoman-Ztiniga et^l2008l ; iPoulton et al.ll2008h . 
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Figure 18. Top panel: Distribution of near-infrared excess source density as a 
function of distance from the center of the Rosette Nebula (NGC 2244). The counts 
are made in sectors of 1 .0 pc in length and counts in each sector have been scaled and 
normalized to the area and counts in the central 1.0 pc circle in NGC 2244. Labels 
indicate the approximate locations of embedded clusters as well as the main 're- 
gions' of the complex, defined by Blitz & Thaddeus (1986). Bottom panel: averaged 
near-infrared excess fractions in each of the cluster groups defined from t he top plot 
appear to increase with distance from the Rosette Nebula. From .Roman-Zuniga et al.l 
(|2008). 
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4. Regions of Particular Interest 



4.1. AFGL 961. Cluster PL06 

Particular attention in the literature has been given to cluster PL06, associat ed with 
the br iaht infrared source AFGL961. The source was originally identified by ICohej 
(Il973h . who also obtained its spectral energy distribution from 2 to 18 /im. This 
distribution was lat er extended to include mea surements at 53, 100 and 175 /xm by 
Harvey et all (Il977h . iBlitz & ThaddeusI (Il980l) found that the CO emission towards 
the source is strongly self absorbed and presents broad wings, evidence of a violent 
interaction of the c entral source with the molecular envel ope. These characteristics 



were confirmed by iLada & Gautier I (Il982h: iLorenl (|l98lh; ISchneider et~al] (Il998al ): 



Heyer. WiUiams. & BruntI (l2006h and iRoman-Zufiigal tOO&l. The powerful CO out 



flow that generates the wings was measured to have a gas kinetic ener gy close to 10^^ 
ergs, and is moving almost 19 M0 of material ( Lada & Gautier |[l982h . The outflow is 
bipolar and spectral signatures revealed that the outflow is in fact caused by an ionized 
wind. The self absorption is possibly caused by a foreground laye r of co lder material, 
and in fact, infrared spectra at 1 /im and 1 mm obtained by ICoxl (Il989i) revealed the 
presence of ice bands towards the source, confirming the effect. 




Figure 19. A FLAMINGOS near infrared image overlaid with contours of 
^^CO(2-l) outflow (line wings) emission towards cluster PL06 which hosts the mas- 
sive PMS binary AFGL961. The solid contours show integrated emission from 5 to 
8 km-s^^, and the dotted contours show integrated emission from 17 to 20 km-s^^. 



The central source in AFGL961 is a B type PMS binary with a separation of ap- 
proximately 0.05 pc, hence the common nomenclature AFGL961a and AFGL961b, 
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which refer to the eastern and western compon ent, respectively. Near infrared observa- 
tions ( Aspinlll998l :lRoman-Zuniga e t al.ll2008 ) revealed that approximately another 30 
young sources are associated with this cluster, and lAspin (1998) was able to demon- 
strate that at least 9 of these sources were associated with shock-like features which 
corroborated the strong interaction with the local medium. 

Figure [l9]combines a FLAMINGOS image with ^^C0(2- 1 ) data from the lRoman-Zuniga 

(I2006h survey in the form of integrated intensity contours. Self absorption is evident 
even in the averaged spectra, and the red and blue wings appear to confirm the exis- 
tence of two different components of the gas outflow, although none of them aie highly 
collimated. 



4.2. Jets, Knots, Elephant Trunks and Globules 

The [SII] survey of lYbarra & PhelpsI (2004) revealed 13 regions of excited gas, 3 as- 
sociated with external excitation from winds of NGC 2244 sources and 10 associated 
with winds from young sources in the embedded clusters of the molecular cloud. Most 
of the main clusters host one of these [SII] features, all reminiscent of HH objects and 
thus suggestive of strong interaction of young sources with their surrounding medium 
or at least of a high frequency of collimated outflows f rom young sources. Counterpart 
near infrared H2 observations (iPhelps & Ybarra|[2005 ) revealed that one source in the 
[SII] study, RMC-C (06315690-1-0419026) could be in fact a Class I object associated 
with a large outflow designated HH 871. The source is very close to cluster PLOl and 



could be associated with it. See Figure 120 



A seri es of studies by iMeabum & WalshI (Il986h . IClayton & MeabumI (1 1995 b and 



Clayton et al. ( 1998) were dedicated to the investigation of a group of dense dense 
nebulous knots located at the south eastern edge of the HII cavity i n the nebula area at 
approx. {a, 6) = (6''29™36^ +5'^00W, J1950). In the first paper. iMeabum & Walsh 
obtained Hq;-i-[NII] filter photographs of the central cavity areas and combined them 
with echelle spectroscopy data. They suggested two possibilities for the origin of the 
knots: a) they are HH objects propulsed by T Tauri sources with bipolar activity, or b) 
they are bow shocks arou nd dense globules ove rtaken by large wind driven shells. 

In the second study, IClayton & MeabumI generated a data cube of [OIII]5007 A 
profiles obtained by scanning the knot areas with a multi-slit echelle spectrometer. They 
noticed that the largest knot, labeled C, is indeed very bright in this line, and suggested 
again that the feature might originate from a low mass young star. In their third study. 



Clayton et all studied again the prominent knot C. The combined observations plus K 



band images identified candidates for the driving sources, and proposed that the top 
group of knots — which appear to coincide with K band sources — are all associated 
with a giant Herbig-Haro flow. The large scale of the corresponding outflow would 
be consistent with the idea of a large jet associated with the source RMC-C (see Fig- 
ure I20I ). thus confirming the existence of parsec scale flows in the Rosette. Several 
FLAMINGOS sources lie in the area and the knots are baiely visible in J and H, so 
a very interesting follow up wo uld be to determine if any other faint sources are also 
associated with the knots. Also, IMeabum et al. 1 (120051) found a possible microjet asso- 
ciated with a low mass young star, and most likely observable due to the Ly-a emission 
of the n ebula. 

O ( 2003 ) and lLi & Rectoil ( 20041) reported the discovery of two a cometary jet sys- 
tems with external irradiation. These two objects were recently catalogued as HH 889 
and HH 890 (see http ://casa. Colorado .ed u/hhcat/, ) . These jets are immersed in the pho- 
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Figure 20. Large outflows and knots in the Rosette Complex environment. The 
left side panel is a composite [SII]-2MASS image of object HH 871 a nd its asso- 
ciated source Rosette Molecular Cloud-C from Phelp s & Ybarral ( 120051) . The right 
panel show s a group of c onspicuous knots in the inner western edge of the HIT re- 
gion. From lClavton et all (,1998,) . 



toionized region and thus are awash by heady UV radiation, which apparently results 
in an interesting set of bow shock structures. Given the age of NGC 2244, there is an 
upper hmit of 1-2 Myr in the timescale of the jet production. Interestingly, the objects 
associated as the driving sources present no infrared excess Li (200j)^ which suggests 
a very rapid disk dissipation, with timescales of 10 ^ to 10^ yr (Li e t al.ll2007 ). leading 
to a rapid transition from CTTS to WTTS stages dLi & R ector 2007|r These aspects 
pose interesting questions about the environmental dependence of the mechanisms of 
jet formation and also suggest that CTTS and WTTS objects can be spatially mixed in 
regions of massive star formation. 

At the northwestern edge of the central cavity in the Ros ette Neb i ila, the re is a large 
group of prominent dusty pillar structures or elephant trunks faerbigl ( i974b . see Figure 
EB. S chneps et al. 01980) observed this region in CO and ^^CO (J=1^0), and found 
that the morphology of the molecular gas emission in the area follows very closely the 
optical outline of the pillar features. A very interesting fact is that the physical prop- 
erties of the globules associated with the larger pillar structures are not different from 
those of globules in clouds not associated with HII regions. Individual globules in the 
area have H2 densities of ~ 10^ cm~^ and temperatures of ~ 10 K. The trunks have 
velocity gradients of up to 1.1 km-cm~^(')^^, oriented along their direction of elon- 
gation and these gradients might be responsible for the morphology. The "stretching" 
of the trunks might have a time scale of ~3 to 6 x 10^ yr, which is consistent with the 
dynamical age of the HII region given an expansion rate of 20 km-s^^. 

Later, the studies of ICarlqvist et al.l (1998; 2002[) were dedicated to study the in- 
ternal structure of the elephant trunks. By constructing an extinction map of the largest 
trunk from a B-band CCD image, they were able to identify a clear helical structure, 
which they suggest might originate from an alignment with a primordial helical mag- 
netic field (magnetic rope model). If true, this magnetic field would be present before 
the HII region was formed, and it could be responsible for maintaining the trunk-system 
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and would explain the observed velocity gradient even after erosion from the winds of 
NGC 2244. 




Figure 21. The northwestern edge of the cavity in the Rosette Nebula showing a 
complex system of pillars, trunks and small globules. Credit: Ignacio de la Cueva 
Torregrosa. 
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Gonzalez-Alfonso & Cemicharol (1994) observed 4 of the densest 'teardrops' lo- 



cated at the tips of the elephant trunks. Using the IRAM 30m telescope they found 
significant emission in 1^00(1-0 and 2-1), ^^COCl-O and 2-1) and CS(2-1 and 3-2). 
They confirmed kinetic temperatures for the globules between 15 and 20 K and masses 
between 0.02 and 0.5 M0. The ^^CO line profiles had significant outflow (wing) com- 
ponents, possibly indicating the pressure ex erted by the ionized gas in the HIT region 
on the neutral gas of the teardrops. Recently. iGahm et all (120071) made a complete sur- 
vey of the entire chain of dark trunks in the northwest edge of the central cavity. They 
observed 23 fields in Ha using the 2.6m Nordic Optical Telescope, obtaining images of 
filaments and globules at a resolution of 0.188 pix^^ (about 300 AU at d = 1600 pc). 
The observations revealed a total of 145 individual small globules, many of them lo- 
cated in isolation near large pillars. They were able to obtain density profiles for the 
globulettes by converting from Aa extinction, and showed that most objects agreed well 
with r^^'^ configurations typical of protostellar conditions. The total densities of the 
globules varied from 3 x 10'^ to 1 x 10^ cm^'^ and values of total mass assuming uniform 
density were found to agree well with those obtained for the few counterparts from the 
survey of Gonzalez- Alfonso & Cemicharol . fCahm et al . found that many objects had 
radii under 10,000 AU, with a distribution peak around 25,000 AU. The distribution of 
mass for the globulettes above a detection limit of 1 M j peaked at about 13 Mj, with 
only a few objects exceeding 60 Mj. No near-infrared counterparts have been found to 
be convincingly associated with the globulettes, possibly indicating that none of these 
objects is a genuine proplyd or a protostar. Virial analysis yields the result that above 
60% of the globulettes have conditions for gravitational contraction. Furthermore, the 
free-fall times were estimated to be of a few times 10^ yr, or about one order of magni- 
tude shorter than the photoevaporation times, estimated to be of the order of 3 to 4 Myr. 
This large photoevaporation timescale is surprisingly large in the harsh environment of 
NGC 2244, but may be due to a process of slow compression that follows the erosion 
of the outer layers by photoevaporation. One of the most interesting conclusions from 
this study is that the small Rosette Nebula globules may be the predecessors to brown 
dwarfs or planetary mass objects, reinforcing the idea that these kind of low mass ob- 
jects can form in isolation following a process of contraction similar to stars, and not 
only by ejection or disruption of protoplanetary disks. 



5. Final Remarks 

After 80 years of studies, the Rosette Complex is known as one of the most important 
astrophysical laboratories for the study of the many aspects of star formation. 

A solid basis for the investigation of the Rosette Complex has been successfully 
set. We have now a very detailed picture of the physical properties of the HIT region 
and the molecular cloud, and we have located the majority of the star s that have formed 
across the complex. The Rosette has been actively forming stars for the last 2 to 3 Myr 
(based on the age of NGC 2244) and it will probably continue doing so, as there is still 
sufficient molecular material (> 1 x 10^ M0), placed in a heavily stimulated environ- 
ment: The large photodissociation region excavated by the OB association NGC 2244 
is in clear interaction with the adjacent molecular cloud, which after several years of 
investigations has revealed an abundant embedded population. Also, it appears that an- 
other population of young, low mass stars could form in the surroundings of the central 
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cavity (NGC 2237, REFLIO and the small globules, if they are truly predecesors to low 
mass objects). Could these be evidence of triggered formation? 

The formation of embedded clusters and associations in the adjacent molecular 
cloud was clearly affected by the interaction with the nebula, but we still have to under- 
stand the extent of that interaction. The different cluster formation episodes appear to 
suggest that the current star forming properties of the Rosette Molecular Cloud could 
be related to the early evolution of the cloud, and not only to the interaction with the 
HII region. 

In the near future it will be necessary to investigate the possible differences in 
the mechanisms of cluster formation, the variations in the star formation efficiency 
across the cloud, the observed properties of clusters and the definitive role of the local 
environment. All of this information will be of enormous value for the global study of 
star formation. 

Now that the second generation of clusters in the Rosette has been identified, it is 
expected to be studied in even greater detail using infrared spectroscopy. Such investi- 
gations will be of great help to constrain the relative ages of the embedded populations 
and to reconstruct the history of formation in the region. Also, the near-future devel- 
opment of observational technologies and continuous improvements in the modeling of 
the interstellar medium and formation of clusters are promising, to say the least, for one 
of the most notorious targets in modem astronomy. 
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